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1 .  SUMMARY 


Part  of  Phase  I  of  the  work  statement  of  Contract  N00014-72-C-0335 
has  been  completed.  This  includesfthe  development  of  a  ray  tracing 
program  for  calculating  the  sound  intensity  from  sources  in  the 
ocean  with  surface  and  bottom  reflections.  Salinity  and  temperature 
gradients  have  been  taken  into  consideration  to  calculate  the  sound 
velocity  and  absorption  loss.  The  sound  intensity  is  calculated  by 
energy  considerations  included  in  the  ray  tube  of  adjacent  rays. 

A  formula  has  also  been  developed  for  calculating  the  intensity  along 
a  ray  without  considering  the  variation  in  the  ray  tube  area.  This 
involves  the  calculation  of  the  radius  of  curvature  of  the  wave  front 
along  the  ray. 

Future  plans  on  ray  tracing  call  for  the  development  of  techniques 
for  calculating  sound  intensities  in  the  caustic  region  and  higher 
order  corrections. 
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2.  DISCUSSION 


ENERGY  ALONG  A  KAY  PATH 

Consider  two-dimensional  ray  propagation  in  the  (x,y)  plane  with  sound 
velocity  c  =  c(x,y) 


y 


X 


Fig.  1 

Let  s  be  the  distance  measured  along  a  ray  and  let  n  denote  the  local 

orthogonal  direction  as  shown.  If  0  is  the  slope  angle  of  the  ray,  the 

•>  ->■ 

(s,  nj  vectors  are  obtained  by  rotating  the  (x,  y)  vectors  counterclock¬ 
wise  through  0. 


The  curvature  1/R  is  given  by 
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It  follows  easily  that 


c  =  cx  sin0  -  cv  cos 0  =  -  ~~~ 
ds  x  /  3n 

so  that 

I  =  In  tSn  c)  (3) 

The  rays  are  orthogonal  to  surfaces  of  constant  phase;  for  brevity,  such 
surfaces  will  be  termed  wave  fronts.  Denote  the  local  curvature  of  a 
wave  front  by  (1/RW);  taken  as  positive  n  the  rays  are  diverging.  The 
immediate  purpose  is  to  compute  d/ds(l/Rw). 

For  convenience,  introduce  an  orthogonal  curvilinear  coordinate  system 
with  §  constant  on  a  ray  and  n  constant  on  a  wave  front. 


Fig.  2 
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Thus, 


d  ,1  ,  ,e§  ..  cosG 

j—  Ct; — J  =  ( —  cosGj  - 

ds  lRw  y§  n  x 


C0S6  [  6k„cos0  -  -  -  0§0nsin0]  (6) 


V§  §n 


§ 


Similarly, 


4-  (i)  =  Soso  .  Msncose 

dn  vIr  xny§  §n  x 


0§Onsin0]  (7) 


combining,  we  obtain 


— r~i  =  ~  r— i  cosS  x§n  y§n> 
ds  Rw  ‘  dn  LRJ  +  R  “  RWJ 


(8) 


Next,  differentiate  the  first  of  Equations  [4)  with  respect  to  §  and 
the  second  with  respect  to  n  j  then  solve  the  resulting  equations  for 
y^^  and  Xg^  to  obtain  (with  the  help  of  Equations  (5}) 

xnT§  _ tan0> 

y§n  ~  cos0  lRw  ~  R  ' 

(9) 

xnT§  ,_1_  tan0, 

§n  cosO  ^~R  "  Rw 


Substitution  into  Equation  (8)  now  gives  the  desired  result; 


•L  (i_)  * 

ds  '■R  J  np  dp 


dnvRJ 


R2  R2 
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(10) 


To  obtain  a  convenient  formula  for  d/dn(l/R),  we  observe  that  Equation  (3) 
yields  the  following. 
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^  =  ;kI(£n  c)x  sine  "  (£n  c)y cos0] 

d0 

=  [ (£n  c)x  cose  +  (£n  c)  sine] 

-  [(£n  c)  .  sin20  -  2(£n  c)  sine  cose  +  (£n  c)  cos2e] 
xx  y  y 

=  ^  [t(£n  c)]  '  S^[£n  c]  (11) 

where  d2/dn2  is  defined  in  the  "straight  line"  sense  via  the  progression 
between  the  last  two  lines  of  Equation  11.  Carrying  out  the  differentiation 
of  (£n  c)  gives  the  alternative  form 


d  fl  •.  _  1  dc  1_  d2c _ 1_ 

ds"  lRwJ  "  cRw  ds  c  dn2  Rw2 


(12) 


where 


4^-  =  c  sin2e  -  2c  „  sine  cosG  +  c  cos2G  (IS) 

dn2  xx  xy  yy 

Suppose  now  that  the  source  is  located  at  some  point  on  the  (negative) 
y  axis.  At  each  point  along  a  ray  we  know  the  range  x,  the  slope  angle 
e,  and  the  wave  front  curvature  1/RW  (as  a  result  of  Equation  (12)). 

Let  F  denote  the  intensity  (energy  rate  per  unit  area)  along  a  ray.  It 
now  follows  from  simple  geometry  that 


I  dF 

"F  ds 


cose  l 
x  R, 


(14) 
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If  A(x,  y)  denotes  the  acoustic  attenuation  in  the  water,  then  the 
final  equation  becomes 


+jL.A(x,y) 

F  ds  x 


(15) 


In  a  program,  it  is  worthwhile  to  record  each  term  in  Equation  (15) 
separately;  thus,  F  is  divided  into  the  three  terms  Fg,  Fjv',  and  F^ 
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via 


F  =  Fr  •  Fjy  *  Fa 

so  that 

1  dFR  _  cosO  1  dx 

-Fr  ds"~  ~  x  x  ds 


1_  _  1_ 
Fw  ds  ~  Rw 


1_  djA 

fa  ds 


A(X,  y) 


(16) 


The  first  of  these  equations  yields 
-  const. 

FR  =  -IT- 

when  the  constant  is  chosen  to  be  the  intensity  of  the  source  at  unit 
distance  (i.e.,  s  =  1),  along  the  chosen  ray;  thus,  F^  and  FA  are  each 
unity  at  that  point. 

A  knowledge  of  F^  has  an  interesting  physical  interpretation.  Consider 
two  rays  emanating  from  the  source  at  angle  0q,  at  an  incremental  angle 
d0Q  apart.  At  unit  distance  from  the  source,  F^  =  1.  At  a  terminal 
point  xt,  where  the  slope  angle  is  and  the  value  of  F^  is  F^,  we  can 

use  the  fact  that  F^  is  inversely  proportional  to  the  normal  spacing  dn 
between  adjacent  rays  to  write 

(l)(l)de0  =  Fwt|dn|  (17) 


so  that  | dn/d0Q |  =  1/F^t 
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